Collapse-driven spatiotemporal dynamics of filament formation 
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The transition from spatial to spatiotemporal dynamics in Kerr-driven beam collapse is modelled 
as the instability of the Townes profile. Coupled axial and conical radiation, temporal splitting and 
X waves appear as the effect of Y-shaped unstable modes, whose growth is experimentally detected. 



From light filaments to Bose Einstein condensates 
(BEC), from plasma instabilities to hydrodynamical or 
optical shocks, many nonlinear wave processes lead to 
self-compression with catastrophic increase in peak in- 
tensity, followed by relaxation to a linear state. Uni- 
versal properties of the compression dynamic were dis- 
covered by analyzing the blowup of self-focusing solu- 
tions to the nonhnear Schrodinger equation (NSE) 
described as a self-similar collapse to a smooth and sym- 
metric ground state Though wave-packets evolve in 
the three-dimensional (3D) space, uneven initial condi- 
tions and compression rates typically lead to symmetry 
breaking, featured by a transient compression towards a 
sub-dimensional eigenstate. In light-beam self-focusing, 
self-similar collapse to the Townes profile (TP) , or ground 
state of the NSE in two spatial dimensions [3| , has been 
recently observed As to the subsequent expansion, 
3D (space-time) X-waves were shown to capture the ap- 
parent stationarity, the pulse-splitting and the conical 
emission (CE) in the filamentation regime [1,0]. Up to 
date, however, the mechanisms that support the transi- 
tion from the quasi-2D, self-similar collapse to the 3D 
relaxation regime remain unrevealed. 

A possible approach to this problem requires retain- 
ing all highly nonlinear effects relevant to this transition 
phase in a suitably "dressed" NSE. Here we show that the 
collapse-driven "morphology" transition from the domi- 
nant 2D to the 3D dynamics can be described in the 
frame of the bare 3D NSE model by the 3D (space-time) 
instability of the 2D TP eigensate. Up to date, transverse 
instability was investigated only for the ID NSE eigen- 
state 0, @!, which however does not support collapse. 
Our analysis directly applies to light filaments generated 
by fs pulses in normally dispersive Kerr media, but can be 
generalized straightforwardly to the collapse in anoma- 
lous dispersion, the spatiotemporal (ST) instability of 
X^^) solitons and the collapse of sub-dimensional localized 
structures in plasmas and BEC. The results capture the 
key steps of the entire filament evolution from the input 
Gaussian to the final X waves, and explain the coupling 
between axial and CE. These were regarded as two inde- 
pendent processes, until recent numerical investigations 
outlined their simultaneous appearance . Here the cou- 
pling emerges as an intrinsic feature of the instability of 



the TP, whose signature is a couple of Y-shaped unsta- 
ble modes. For the first time, a theory is provided that 
links filamentation, CE, pulse-splitting and supercontin- 
uum generation to the singular behavior at collapse. 
Consider the cubic, dimensionless NSE 

^cA^'-A^,,A-'-^lA + ^\A\'A, (1) 

with A^ jj = for a function A that depends on the 

spatial coordinates (^, rj) and the temporal coordinate r. 
The ground state of the NSE in two spatial dimensions 
is expressed hy A = an(p)exp{ia(^) [p = (^^ -|- ??^)^^^], 
where ao(p) is the TP 0] [Fig. [iKa)], and a ~ 0.2055 for 
the normalization ao(0) = 1. 

Written in the variables z = (no/A:on2/)C, ix,y) = 
{no/kln2lY'^{^,7^), t-k'oz = (nofc(,'/fcon2/)i/2r, A = 
1^1^ A, Eq. C]) is a standard model for the propagation of 
an optical wave packet E — A{x, y, z, t) exp(— iwoi+ifcoz) 
subject to diffraction, normal group velocity dispersion 
(GVD) (/cq > 0) and self-focusing nonlinearity (n2 > 
0) leading to collapse. The ground state describes a 
monochromatic light beam with z-independent trans- 
verse amplitude profile. In the above relations, luq is 
the carrier frequency of the wave packet, / is a charac- 
teristic intensity, n2 the nonlinear refraction index of the 
medium, and fc(a;) — n(uj)uj / c the propagation constant, 
with n(uj) the refractive index and c the speed of light 
in vacuum. Prime signs stand for differentiation with 
respect to lo, and subscripts for evaluation at ljq- 

Following a standard MI analysis, the perturbed TP 

A = {ao(p) + a[u{p)e-'^^+""^ + v*{p)e'^^^"'*^]}e'°'^ , 

(2) 

with (T <C 1, is considered. The ST perturbation will grow 
exponentially with C if k/ < (k/ = Imre). Opposite 
frequency shifts and — f2 characterize the components 
u and V* . For simplicity, only cylindrically symmetric 
perturbations are considered. Introducing ([2]) into ([T]), 
and keeping terms up to the first-order in cr, one gets the 
differential eigenvalue problem 

where H = /dp^ + {1/ p)d/dp + VL^]-a + 2al{p), with 
boundary conditions du{Q)/dp = dv{0)/dp = 0. Given a 
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-[ri^ — 2a + 2k]v are obtained. Their bounded 



FIG. 1: (a) Normalized TP. (b) Real and imaginary parts of 
the relevant complex eigenvalue. 



perturbation frequency fl, if the problem admits at least 
one nontrivial solution (m, v) with an eigenvalue k with 
negative imaginary part, then the TP is unstable un- 
der perturbations at that frequency. The problem ([3|) 
has to be solved numerically. To simplify the procedure, 
note first that if {k,u,v) is a solution of ^ at fi, then 
{k*, u*, V*), (— K, V, u) and (— k*, v* ,u*) are solutions too, 
whose eigenvalues are all reflections of k about the real 
and imaginary axes in the complex K-plane. Then, if k, 
is complex, only two of these solutions represent unsta- 
ble perturbations. Following the convention that kji > 
[kr = Ken) and k/ < 0, these two solutions are (k, u, v) 
and {—K*,v*, u*), that yield the two independent, unsta- 
ble and physically distinguishable perturbations 

since u and v* are generally different and oppositely fre- 
quency shifted in and p-K*- Also, any solution of ([3]) 
with f2 is also solution for — fi. In particular, the unstable 
solutions {k,u,v) and {—k*,v*,u*) for the perturbation 
frequency — O are seen from ([2]) to represent the same 
but permuted perturbations p-K.* and p^- In short, we 
can restrict the analysis of ([3]) to > 0, and count only 
eigenvalues with kr > 0. For each of these with ki < 0, 
two independent, unstable perturbations ^ may grow. 

The problem ^ has been solved numerically for each 
il > after discretization of the differential operators 
on a radial (p) grid of finite size much larger than the 
Townes range d ^ 1.5. Increasing grid points and size 
allowed to control the accuracy of the results. For SI = 0, 
no complex eigenvalue is found, as expected from the 
marginal instability of the TP under pure spatial pertur- 
bations 0]. For each > 0, only one, isolated eigenvalue 
K with K/ < appears, whose real and imaginary parts 
are shown in Fig. [l]Jb). The ground state ao{p) exp{iap) 
is then retrieved to be modulationally unstable under ST 
perturbations The gain —kj is limited to < 1.5 

but without an abrupt cut-off. 

For the eigenmodes u and v associated with the com- 
plex eigenvalue k, some relevant features can be in- 
ferred without resorting to numerical calculation. Ne- 
glecting terms with ag in ([3]) and in H for p —> oo, 
the decoupled equations A^^^u = — [f2^ — 2a — 2k]u and 



solutions for p CX3 are u oc [{Qu + i^u)p] and 
V cx Hq^^ [{Qv +i^v)p], where H^^'' is the Hankel func- 
tion of first class and zero order, and the real quantities 
Qu,v and r„ are defined by 



Qu+iTu=\/^^-2a-2K, Q^+iT„^y/n^ -2a+2K, (5) 

with the convention of taking square roots such that 
r„,„ > 0. Using that H^^\s) ~ V(2/7rs)exp [i{s - 7r/4)] 
for large |s|, ignoring constant factors and algebraical 
decay, the expressions u{p) ~ exp (— r„p -I- iQuP) and 
v{p) ^ exp (— r„p -I- iQvp) are found to describe the dom- 
inant behavior of the unstable eigenmodes for large p. 
Thus, if their decay rates are small {Tu,v ^ 1)7 u and 
V are expected to feature radial oscillations of transverse 
wave numbers Qu,v beyond the TP range. 

The physical quantities featuring the collapsing dy- 
namics are the growing perturbations p^ and p-^* com- 
posed of u and v. In particular, the ST spectrum Q-Q 
of the unstable perturbations at a certain propagation 
distance C will be generally expressed by a superposition 
of the type af,{fl)pi^{Q,n) -I- (7_„. (ri)p_^* (Q, O), where 
Pii{Q,fl) and p-^* {Q,i^) are the spatial Fourier trans- 
forms ofpK{p,n) andp_K*(p, fi), and ^^(ri) and (T_K*(ri) 
are the seeds of each unstable perturbation at each fre- 
quency. In experiments, measurement of the ST spec- 
trum Q-fl (in practice, the angularly-resolved spectrum, 
displaying angles and wavelengths) of a light filament is 
a powerful diagnostic in which axial and CE can be visu- 



alized at once . For the TP, the structure of the Q-fl 
spectrum of instability can be inferred from the above 
asymptotic analysis. The spectrum p^iQ, ^) is composed 
of u{Q, fl) at positive frequency shift and expectedly 
peaked at Qu, and of v*{Q, f^) at negative frequency shift 
—fl and peaked at (if the respective decay rates r„^„ 
are not large). Plotting Qu and Qv given by ^ versus 
frequency in their respective ranges, we obtain the black, 
thick solid curve of Fig. Ufa). This curve is expected 
to locate the regions in the Q-U plane where p^iQ,^) 
takes its highest values. Fig. [Ifa) shows also the values 
of r^, and r„ (black, thin solid curve). Since exponen- 
tial localization of v is weak (F^, <C 1), the branch Q^, 
fitting approximately the linear relation Qy ~ —y/2il, is 
expected to define the locus of maxima of v{Q,^), and 
hence Pk{Q, ^) at down- shifted frequencies to describe an 
actual off-axis, CE. Instead, localization of u is similar to 
that of the TP (r„ ~ 1) and Qu close to zero, which sug- 
gests to associate Pk((5, ^) at up-shifted frequencies with 
an axial emission. Similar analysis holds for p_^*((5,S7), 
CE being now associated with up-shifted frequencies, and 
axial emission with down-shifted frequencies [see Fig. 
[21[a), black, dashed curves]. This characterization of the 
instability spectrum is confirmed by numerical evaluation 
of u and v. Fig. [H^b) shows Pk{Q, ri) with the expected 
features. Its refiection about SI = yields j5_„* (Q, S7). 
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The instabiUty of the TP then reveals the link between 
axial and CE. In this view, the basic ingredients of the 
ST spectrum of a light filament are not the axial contin- 
uum and the X-shaped CE, but two Y-shaped objects, 
each one coupling off-axis emission to axial emission on 
opposite frequency bands. 

Let us compare these results with those for the ST in- 
stability of the plane wave solution A = exp{i() to the 
NSE [3- The problem Q holds for this case if the 
amplitude ao and the nonlinear wave vector shift a are 
replaced by unity, and presents the same symmetries. 
The eigenmodes u and v associated with the most unsta- 
ble perturbations and p-^* at each frequency fl are 
two plane waves with identical transverse wave numbers 
Qu = Qv = Vil^ + 2. The gain — k/ = 1 is indepen- 
dent of fl and hence unbounded. The spectra pi^{Q,Q,) 
and p-K*{Q,^) are then characterized by two identical 
hyperbolas, depicted in Fig. ^a.) as a single green curve 
(see Ref. [12|]), and usually referred to as an X-like spec- 
trum. The slope of the X arms is unity, while that of the 

Y arm is \/2 ( y^fcofc^ against •\/2fcofco physical units) . 
Degeneracy breaking from single X to double Y and 

gain limitation originate from the transverse localization 
of the TP. To see this more intuitively, we interpret insta- 
bility as originated from a Kerr-driven, degenerate four- 
wave-mixing (FWM) interaction in which two intense, 
identical pump waves of frequency = propagating 
along the C direction, amplify two weak, non-coUinear 
plane waves, u and v, at frequencies ±17. Each pump 
wave propagates with an axial wave number shift fciML 
(with respect to a plane wave of frequency SI = 0) due 
to self-phase modulation (SPM). The axial wave num- 
bers of the weak plane waves are shifted by — (5,^/2 
and fl'^ /2 — Ql/2 [as obtained from the hnearized NSE 
([T])] due to material dispersion and tilt, and by fcNL,ti and 
/cnl.d due to cross-phase modulation (XPM). Axial phase 
matching among the interacting waves then reads as 

flV2-Q2/2 + fcNL,« + r!V2-Q'„/2 + fcNL,. = 2fcNL. (6) 

For plane wave pumps, SPM and XPM wave number 
shifts are /cnl = 1, ^nl,« = fcNL.u = 2fcNL = 2 [ll]. 
Axial and transverse phase-matching {Q^ — Ql) must 
be strictly enforced for efficient amplification of u and v, 
which lead to Qu.v — vfP + 2, i.e., to the X-like spec- 
trum of instability of the plane wave. 

Suppose, instead, that the pumps have TPs (/cnl = a)- 
To account for transverse localization, a) we will assume 
that plane waves coUinear with the pump are preferen- 
tially amplified, and b) we take into account that effi- 
cient amplification is also possible with small transverse 
phase- mismatch such that \Qu — Qv\ ^ 7r/d If for 
instance, the plane wave u at +fl is taken as collinear, 
then Qu — 0, and fcNL,ii — 2/(;nl = 2a. Axial phase 
matching ([6]) then requires Q„ 9^ 0, i.e., the plane wave 

V at —fl cannot be collinear. It is then reasonable to 
take its XPM wave number shift as fcNL.u — 0, since the 
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FIG. 2: (a) Thick black curves: Characterization of the ST 
spectrum of the unstable perturbations Pn and p-K* by means 
of the transverse frequencies Qu,v Thin black curves: Expo- 
nential radial decay rates r„,„. Green curve: maximum gain 
curve of the plane wave. Red and blue curves: double-X spec- 
trum, (b) Modulus of the ST spectrum pii{Q,Q,) of p^. The 
amplitudes for different Q are arbitrarily chosen so that the 
energy Ip^fQdQ is independent of Q. 



beam does not remain into the localized interaction area 
With these assumptions, ([6]) leads to the linear re- 
lation Qy ~ — V217, as for the Y-tail in the instability 
spectrum of the TP. Limitation of the gain band finds 
also explanation from the maximum allowed transverse 
mismatch. If Qu — 0, then Qy < ir/d. Taking d ~ 1.5 for 
the TP, we obtain Qy ^ 2, which corresponds roughly to 

^ 1.5, as obtained above. Obviously, if the collinear 
plane wave u is taken at —fl, similar considerations lead 
to the reflected Y wave. 

It follows from our analysis that the onset of down- 
shifted axial emission must be accompanied by upshifted 
CE, and vice versa, and that these two events may occur 
independently. Our simulations and experiments support 
this interpretation. Figure [3U a) shows the ST spectrum 
of a TP that is not a CW but a pulse. This "perturbed" 
TP is moreover asymmetric in time by choosing different 
durations for the leading and traifing parts (Ar = 0.87 
and 29, respectively). In Fig. ^h) the difference between 
the propagated spectrum under the NSE (1) (distance 

= 7.6) and the input spectrum is shown in order to 
visualize the newly generated frequencies. A single Y- 
shaped structure is formed, with a slope of the conical 
part fitting to \/2 (thick line). 

In the experiments, we used a 15 cm- long fused sil- 
ica sample as Kerr media. 200-fs-long pulses at 527 nm 
delivered from a 10 Hz Ndiglass mode-locked and regen- 
eratively amplified system (Twinkle, Light Conversion), 
were spatially filtered and focused with a 50 cm focal 
length lens. The pulses then entered into the sample, 
whose input facet was placed at 52 cm from the lens, 
and formed a single filament for input energies E ^ 2 ^J. 
Angularly resolved spectra of the filament at the output 
facet and from single laser shots were measured with an 
imaging spectrometer and a CCD camera, as described 
in [ll|. At 2 fxJ in fused silica the filament is formed 
just before the output facet of the sample, the Y-shaped 
(blue axial, red conical) spectrum of Fig. being then 
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FIG. 3: ST spectrum (a) of input asymmetric pulsed TP (see 
text), (b) after propagation a short distance (subtracted of 
the input spectrum). Angularly resolved spectra of filaments 
in fused silica at {c) E = 2 /ij and (d) E — 3 fiJ. 



observed. At 3 fiJ, the filament is formed closer to the 
input facet. The double Y spectrum of Fig. ^d) then 
corresponds to a longer filament path within the sample. 
The faster growth of one of the two unstable perturba- 
tions at the initial stage of filamentation supposes some 
unbalancing in their seeds. In fact, an expected insta- 
bility seeding mechanism, as self-phase modulation, gen- 
erates preferentially axial, blue-shifted frequencies in the 
presence of plasma and third-order dispersion. 

ST instability of the TP also incorporates a mecha- 
nism of temporal splitting. If the perturbation p„ (e.g.) 
is seeded coherently at different frequencies, its growth 
leads to the formation of a pulse, whose spectrum is 
increasingly peaked at the maximum gain frequencies 
Om — 0.365 for u and —flivi — —0.365 for v. For u 
and V with respective axial wave number shifts kr + a 
and —kr + a, the inverse group velocity (in the 

frame moving with the group velocity of a plane pulse at 
the TP frequency) will be [dKu/dn\fij^,j] ~ 0.5 for u, and 
[d(— /«fl;)/(ir2|_n„] — 0.5 for v. Being equal, the Y-wave 
propagates as a whole with a well-defined group veloc- 
ity. For the perturbation, or reflected Y-wave, the 
group velocity is the opposite, and the group mismatch 
between the two Y- waves is — l/'^^p^''^* — li which 

in physical variables yields, •\/fcofcQr7,2//rto, in agreement 
with the observed dependence in filamentation on pump 
intensity and material properties [6]. 

If the pump is not strictly monochromatic but a long 
pulse, the Y- waves will leave it at a certain stage, ceas- 
ing then to grow. The Y-wave behaves then as a 
linear wave of zero mean frequency and zero axial wave 
number shifts, since u and v, with opposite frequency 
shifts ±S7m, have also opposite axial wave number shifts 
f^M/2-OS/2 nlj/2 and nlj/2-Ql/2 ~ ~nlj/2. The 
peculiarity of the linear Y-wave is that the u part expe- 
riences the normal GVD +ri^/2 of the medium, while 



the V part experiences a net anomalous GVD —57^/2 
as a result of material and angular dispersion. This 
makes the u and v parts of p„ to continue to depart from 
the pump at identical group velocity [(i(^2^/2)/(i^2|oJ^,] = 
[d{-n'^/2)/dn\-n^,] = ^ 0.365 (again, the group 
velocity of the p-^* Y-wave is the opposite). 

Extending our analysis, we may venture an explana- 
tion to the fact that two X- waves are commonly observed 
at later stages of propagation Consider, within 

the FWM approach, the possible effects on the instabil- 
ity spectrum of the strong temporal localization of the 
pump, as may take place upon (possibly multiple) split- 
ting. For a spatially and temporally localized pump of 
central frequency £1 = 0, axial wave vector shift /cnl due 
to SPM, and propagating at a group velocity v^-^^ dif- 
ferent from that of a plane pulse at = (as for an 
Y-wave), new plane waves u and v at opposite frequen- 
cies ±r2 are expected to be preferentially amplified if in 
addition to axial phase matching [Eq. ([6])], the veloc- 
ity of the group formed by u and v matches the velocity 
of the pump. Equating the inverse beat group velocity 
[in^/2 - Ql/2) - {n^/2 - Ql/2)]/2n between u and v to 
l/v'^^\ we obtain _ q2 _ /^Q/v^s\ Phase and group 
matching then yield Qu,v — ^/il^ — {±n)/v^3) — 2fcNL 
{+n for u, —il for v), which is the so-called dispersion 
curve of a frequency-gap X-wave mode 0, [l^. If this 
process takes place for the two Y-waves with group ve- 
locities — ±Qm, the two X-waves of Fig. (^^a, red 
and blue curves) with opposite frequency gaps are gener- 
ated from the two Y-waves. On nonlinearity relaxation 
(^NL 0) at increasing distances, one branch of each 
X-wave is seen to pass through the point Q = = of 
the spectrum, as frequently observed fs^, '6] . 

In conclusion, Y-waves feature the MI of the TP and 
constitute the missing links between self-similar spatial 
self- focusing and ST collapse driven dynamics. Y-waves 
model the coupling between axial and CE, and allow us 
to interpret temporal splitting, X-wave formation and fi- 
nal relaxation into linear X-waves. These results provide 
a unified view of ultrashort pulse filamentation and are 
relevant to all systems involving ST coupling and nonlin- 
earity such as x''^'' solitons, BEG, etc. 
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